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Abstract 

The calculations of nuclear matrix elements of Oi^/5/5-decay is a 
challenge for nuclear physics. We are discussing here a model in- 
dependent method, which could allow to test the calculations. The 
method is based on the factorization property of the nuclear matrix 
elements and requires observation of neutrinoless double /3-decay of 
several nuclei. 



1 Introduction 

The discovery of neutrino oscillations in the atmospheric Super Kamiokande 
, solar SNO j2] , reactor KamL AND j3] , accelerator K2K [4 and other neu- 
trino oscillation experiments 012113111 is one of the major recent discoveries 
in the elementary particle physics. All neutrino oscillation data with the 
exception of the data of the short baseline LSND experiment are well 
described by the three-neutrino mixing ^ 

3 

l^lL = ^ UiiUiL . (1) 
1=1 

Here U is PMNS IT^ mixing matrix, is the field of neutrino with mass 
nii and {I = e, fi, r) is neutrino field which enter into the standard charged 
and neutral currents (so called flavor neutrino field). 

"'^A report at the hiternational Workshop "Neutrino Telescopes", Venice, February 22- 
25, 2005 

^The LSND result will be checked in the near future by the running MiniBooNE ex- 
periment jlOj . 
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From the analysis of the SK atmospheric neutrino data for oscillation 
parameters |Am32| and sin^ 2^23 the following 90% CL ranges were found 

1.9 < |Am^2| < 3 ■ 10"^ eV^; sin^ 26*23 < 0.9 (2) 

From analysis of the data of the solar and the reactor KamLAND experiments 
for oscillation parameters Arn^i and tan^ 612 it was found 

Ami, = (8^^) ■ 10-^ eV^; tan^ 9^ = OAbt'oZ (3) 

From the exclusion plot obtained from the analysis of the data of the re- 
actor CHOOZ experiment ^2] for the parameter sin^ ^13 the following upper 
bound can be inferred 

sin^ ^13 < 5 • 10^2 (90% CL). (4) 

An information about the absolute values of neutrino masses can be ob- 
tained from experiments on the measurement of the high-energy part of the 
/9-spectrum of tritium and from cosmological data. From the data of the 
Troitsk ^3] and Mainz tritium experiments it was found 

mo < 2.2 eV, (5) 

where mo is the mass of the lightest neutrino. From analysis of the cosmo- 
logical data for the sum of the neutrino masses the upper bound in the range 

^m^ < (1-2) eV. (6) 

i 

was obtained (see [IS])- 

The nature of the massive neutrinos (Dirac or Majorana?) is at present 
unknown. The progress in the understanding of the origin of the neutrino 
masses and mixing strongly depends on the answer to this fundamental ques- 
tion. In particular, the proof that massive neutrinos are Majorana particles 
would provide a strong argument in favor of the famous see-saw mechanism 
which, apparently, is the most natural mechanism of the generation of small 
neutrino masses. 

■'The neutrino mass-squared difference is determined as follows: Am.^^. = mf — m^. 
In the case of normal mass spectrum neutrino masses are labeled in such a way that 
TOi < 7712 < iTT-s- In this case A777|2 > 0. In the case of inverted neutrino mass spectrum 
neutrino masses are labeled as follows 7713 < 77ii < m2. In this case Am^2 < 0- 
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2 Neutrinoless double /3-decay 



To reveal the nature of the massive neutrinos it is necessary to study processes 
in which the total lepton number L = Le + + is not conserved. The 
search for neutrinoless double /3-decay (Oz//3/5-decay) of some even-even nuclei 

{A, Z) -^{A,Z + 2) + e- + e~ (7) 

is the most sensitive method of the investigation of the nature of the massive 
neutrinos Vi. 

If Vi[x) satisfies the Majorana condition 

v,{x)=v'i{x)=CvJ{x), (8) 

(C is the matrix of the charge conjugation) the process ((Tj) is allowed. In 
this case neutrinoless double /3-decay is a process of the second order in the 
Fermi constant Gp with virtual neutrinos Vi. The half-life of the process is 
given by the following general expression (see reviews [T7] ) 

^ A = \^ee? \M{A, GiEo, Z). (9) 

Here 

m,, = J2ulm, (10) 

i 

is the effective Majorana mass, G{Eq, Z) is known phase-space factor {Eq is 
the energy release) and M{A, Z) is nuclear matrix element (NME). For light 
neutrino masses NME do not depend on rrii. 

The results of many experiments on the search for OujSf] -decay are avail- 
able at present (see flBj). The most stringent lower bounds on the half-life of 
the Ouf] l3- decay were obtained in the Heidelberg- Moscow ^^Ge experiment 
[T^ and in the recent CUORICINO ^^°Te cryogenic experiment pUj : 

Ti/2(™Ge) > 1.9 ■ 10^^ y (Heidelberg - Moscow) 

Ti/2(^^°Te) > 1.8 • lO^V (CUORICINO). (11) 

Taking into account different calculations of nuclear matrix elements, for the 
effective Majorana mass \mee\ from these results the following upper bounds 
can be inferred 

\mee\ < (0.3 — 1.2) eV (Heidelberg — Moscow) 

keel < (0.2- 1.1) eV (CUORICINO) (12) 
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Many new experiments on the search for the neutrinoless double /3-decay 
of different nuclei are in preparation at present (see |2I])- In these future 
experiments two order of the magnitude improvement in the sensitivity to 
|mee| is expected : 

\mee\ — a. few -10"^. (13) 

3 Effective Majorana mass 

The effective Majorana mass niee is determined by neutrino masses rrii and 
elements Uei. In the standard parametrization we have 

Uei = cos6'i3 cos 6^"^ Ue2 = COS 6'i3 slu 6*"^ U^a = siu 6'i3 e*"^ (14) 

where ctj are Majorana CP phases. 

From neutrino oscillation data we know the values of neutrino mass- 
squared differences jArngg] and Am^i and the angle 623 (see (0) and 0). 
We know also that the angle ^13 is small (see (0}). We do not know the 
character of neutrino mass spectrum (normal or inverted), the mass of the 
lightest neutrino mo, which determine the absolute values of neutrino masses, 
and Majorana phases. 

Let us consider three standard neutrino mass spectra (see [221 ) 

1. Hierarchy of neutrino masses 

■nil <^ 1712 <^ w-s- (15) 

In this case we have 



sin^ 012 \/ Amii + sin^ 6*13 \/ Aml^ e^'"^^ 



< ( sin^ 9i2 \l Am|i + sin^ 6113 \l Am|2 i , 



(16) 



where 032 = a?j — a2. Using the values (jS)) and of neutrino oscillation 
parameters from (fTH|) for effective Majorana mass we find the upper 
bound 

|mee| < 6.4 ■ 10-3 (17) 
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Thus, in the case of the neutrino mass hierarchy the predicted upper 
bound of \mee\ is significantly smaller than the sensitivity of the most 
ambitious future experiments. 

2. Inverted hierarchy of neutrino masses 

^ mi < 1712- (18) 

In this case we have 



|mee| ~ y |Am|i| (1 - sin^2 6'i2 sin^a2i)2 (19) 

From this relation we find 

■\J |Am|J cos 2 < \mee\ < \j \^m\i\ (20) 

Let us notice that the bounds in this inequality correspond to the case 
of the CP conservation in the lepton sector: the upper (lower) bound 
corresponds to the case of equal (opposite) CP-parities of P2 and Vi. 

Thus, in the case of the inverted neutrino mass hierarchy from neutrino 
oscillation data we can predict the upper and lower bounds of the pos- 
sible values of the effective Majorana mass. From (j2)), and pO|l we 
obtain the range 

1.0 ■ 10"^ < |mee| < 5.5 ■ 10"^ eV (21) 

The values of \mee\ in this range apparently will be reached in future 
Ou(3P- experiments. From ()19|1 for the parameter sin^a2i, which char- 
acterizes CP violation in the case of the Majorana neutrino mixing, we 
have ^ 

sin^ a2i = — I 1 - , I (22) 

sin2 2^i2 V l^^lilJ 

Hence, the observation of the neutrinoless double /3-decay with \mee\ 
in the range (pT|l could allow to obtain an information about Majorana 
CP phase difference. 

3. Quasi-degenerate neutrino mass spectrum 



mi < 7722 < "^3 (?^3 < rrii < 1712); mo ^ -v/|Am|2| (23) 
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Neglecting small contribution of sin^ ^13 for the effective Majorana mass 
we have in this case 

\mee\—mo (1 — sin^ 2 sin^ 021)^ (24) 

In the future tritium experiment KATRIN [SB] the sensitivity mo — 0.2 
eV is planned to be reached. If the mass of the lightest neutrino mo 
will be measured in this experiment for the effective Majorana mass 
from and we obtain the range 

0.23 mo < \mee\ < mo- (25) 

On the other side if neutrinoless double /3-decay will be observed in 
future experiments with the value of \mee\ which is larger than the 
upper bound (j21|) for the mass of the lightest neutrino we will have the 
bound 

\mee\ <mo< 4.4 \mee\- (26) 



4 Nuclear matrix elements 

The observation of neutrinoless double /3-decay would be a direct proof that 
Ui are Majorana particles. As we have seen, the determination of the effec- 
tive Majorana mass would allow to obtain an important information about 
character of the neutrino mass spectrum, mass of the lightest neutrino and, 
possibly, Majorana CP phase difference. 

However, from the measurement of the half-life of Oz^/3/?-decay only the 
product of the effective Majorana mass and nuclear matrix element can be 
determined. The calculation of NME is a complicated nuclear problem (see 
|.17j). Two models are commonly used: Nuclear Shell Model (NSM) and 
Quasiparticle Random Phase Approximation (QRPA) with numerous modi- 
fications. These two approaches are based on different physical assumptions. 
As a result different calculations of the same NME differ by factor 2-3 or 
even more . 

We will discuss now a possible method of a model independent test of 
NME calculations j21I. We will use only the general factorization property of 
matrix elements of the OiyjSP decay. Namely the fact that the matrix element 
of the process is a product of the effective Majorana neutrino mass, which 
is determined by neutrino masses and mixing, and nuclear matrix element. 
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which (for hght neutrinos) does not depend on neutrino masses. From (jU)) we 
find the following relations between half-lives of the Oupf]- decay of different 
nuclei: 

Ti/2(Ai, Zi) = X{Au Zi; A2, Z2) ri/2(^2, Z2) = 

Z,; A3, Z,) T,/2{A,, Z3) = .... (27) 

Here 

X(A 7- A .. \M{A,,Z,)\'G{E,,Z,) 

x{A, z., A„ z,) - j^i{a-zWg{e^ ^''^ 

The coefficients in relations (j^7|) have to be calculated. If Oz//3/?-decay of 
different nuclei will be observed in future experiments and relations ()27p 
with coefficients X{Ai, Zf, Ak, Zk) calculated in some model M are satisfied 
than the model M is compatible with data (it is obvious that if relations (j7fj) 
are not satisfied the corresponding model must be rejected). This does not 
mean, however, that the model M allows us to obtain the correct value of 
\mee\ from experimental data. In fact if nuclear matrix elements, calculated 
in the framework of different models and are proportional 

\MmM,Z)\=P\MmM^Z)\, (29) 

{l3 is a coefficient which does not depend on {A, Z)) and relations are 
satisfied for the model Ma than obviously they are satisfied also for the model 
Mfc. The values of the effective Majorana mass, which can be determined from 
experimental data in the framework of these two models, are connected by 
the relation 

1,1,1 / N 

\mee\Ma = ^ ImeelM^ (30) 

and could be quite different. 

For the purpose of illustration we will calculate the coefficients in Eq. fl27p 
in three different recent models of the calculation of nuclear matrix elements 
of Oi^/SP-decay: RFSV CS [26] and NSM In the paper [2S], based 
on QRPA and renormalized QRPA approaches, the values of the parameter 
of the particle-particle interaction gpp, determined from the measured half- 
lives of the 2z//3/3-decay of corresponding nuclei, were used. In the QRPA 
calculation [2^] the values of the parameters were determined from the data 
on the /3-decay of the nuclei of the interest for double /?-decay transitions. In 
the paper ^27, the results of the latest NSM calculations were given. We will 
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consider four different nucleus: Ge, ^^''Mo, ^^"Te and ^^^Xe. In the Table 
I we have presented the values of coefficients X{Ai, Zf, A^, Z^) in the case if 
Oz//3/3-decay of ™Ge and one of the other nucleus is observed. 

Table I 

The values of the coefficient X{Ai, Zf, Ak, Zk) obtained with NME 
calculated in |25| (RFSV), in |21 (CS) and in ^ (NSM). 







RFSV 


CS 


NSM 


X 


[i^UTe;^« Ge) 


0.38 


0.13 


0.24 


X 


'i36Xe;76 Ge) 


0.80 


0.07 


0.56 


x{ 


looMo;^^ Ge) 


0.59 


0.17 





If it would occur that the relation (fTTj) with the coefficient, calculated in one 
of the model considered, is satisfied, in this case, as we can see from the Table 
I, other models apparently can be excluded (if accuracy of experimental data 
are better than ~ 30 %) 

Let us stress that this conclusion depends on nuclei for which neutrinoless 
double /?-decay is observed. For example, if Oi^/3/3-decay of ^^°Te and ^°°Mo 
will be observed, in this case we have 

X(i°°Mo;^6 Ge) = 1.5 (RFSV); 1.3 (CS) 

The difference between the values of these coefficients is only ~ 10 %.. Thus, 
if the relation is satisfied for NME calculated in, say, RFSV model |2S] , 
it will be difficult to exclude the model CS [2^1 • The values of the effective 
Majorana mass which can be obtained from the experimental data in the 
case of these two models are, however, quite different: 

|"iee|RFSV - 2.5 ■ |mee|cS- (31) 

The observation of Oi^/?/3-decay of three (or more) nuclei would be an impor- 
tant tool in a model independent approach to the determination of the value 
of the effective Majorana mass |mee| which we have discussed here. 
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5 Conclusion 



The establishment of the nature of the massive neutrinos Pi would have a pro- 
found importance for the understanding of the origin of small neutrino masses 
and neutrino mixing. The investigation of the neutrinoless double /3-decay is 
the most sensitive method which could allow us to reveal the Majorana na- 
ture of the massive neutrinos. Today's limit on the effective Majorana mass, 
which can be inferred from the study of this process, is \mee\ < (0.3 — 1.2) eV. 
In several experiments now in preparation the sensitivity \mee\ — a few 10~^ 
eV is planned to be reached. If \mee\ is measured, the pattern of the neu- 
trino mass spectrum and, possibly, Majorana CP phase could be revealed. 
For that not only Oz//?/5-decay must be observed but also nuclear matrix el- 
ements must be known. Observation of 0z^/?/3-decay of several nuclei could 
provide a model independent method of testing of NME calculations. 
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